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Abstract 



Motivated by a proposal of Daykin [T], we study the wave that propagates along an 
infinite chain of dominoes and find the limiting speed of the wave in an extreme case. 
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> ' 1 Introduction 

00 : 

^ ■ Everyone is familiar with dominoes and has used them for fun. A common game is to arrange 

ps) ■ the dominoes in a row and give a push to the first. This generates a pleasing wave of falling 

. dominoes. The propagation happens at some speed v (not necessarily constant). A qualitative 

O ' discussion for general audiences is given by Walker in [2]. Given the game's simplicity it is 

^ ; perhaps surprising to discover that an exact computation of the speed v is quite difficult. 

' Daykin realized this in the following 1971 proposal [1] to the readers of the SIAM Review: 
• ^ . 

' "How fast do dominoes fall?" 

^ . The "domino theory of Southeast Asia" says that if Vietnam falls, then Laos falls, then 

■ ■ ■ ' Cambodia falls, and so on. Hearing a discussion of the theory led me to wonder about 

the proposed physical problem. The reader is invited to set his or her own "reasonable" 
simplifying assumptions, such as perfectly elastic dominoes, constant coefficient of friction 
between dominoes and the table, and initial configuration with the dominoes equally spaced 
in a straight line, and so on. 

In 1983, McLachlan et al. [3] found a scaling law for the speed v in the limiting case of 
dominoes with zero thickness equally spaced in a straight line. With these assumptions the 
authors found the functional relation 

M""^- (1) 
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Here i is the height of the dominoes, d the spacing between dominoes, and G{x) an unde- 
termined function of x. This relation followed from dimensional analysis of the proble 
McLachlan et al. proceeded to test the formula experimentally using dominoes of heights 
h = AAA5cm and h = 8.890cm. More recently. Banks presented a simplified description of the 
effect His analysis, among other assumptions, assumes a uniform propagation speed and 
conservation of momentum. A different direction was taken by Shaw in a short paper [3] de- 
scribing how to model the domino effect as a computer simulation and use it as an experiment 
in the undergraduate physics lab. 

Equation ([1]) is not an entirely satisfactory solution, of course, containing as it does an 
unknown function. This recognition became the motivation for the present article. In this 
work we develop an expression for the speed v. The result can be cast as a particular forms of 
the scaling function G arising from a particular set of assumptions. 

In order to highlight the basic physics behind the problem, we replace the dominoes by 
massless rods topped with point masses m, as seen in figure [2l A similar analysis in the case of 
dominoes with the shape of parallelepipeds is straightforward, although there are some minor 
differences. 
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Figure 1: A uniform chain of dominoes. 
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Figure 2: A chain of massless rods carrying masses m on top. The masses are indicated as 
finite spheres only for the sake of visualization. 



^This paper should allow the reader to write down a complete list of assumptions needed to reach this 
conclusion. 
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2 The Model 

2.1 The Assumptions 

We shall assume that: 

1. The chain of rods is uniform. This means that all rods are identical and are equally 
spaced along a straight line. Let i be the length of a rod and m the mass on top. 

2. The collisions are head-on. This means that, seen from above, all the rods are and 
remain aligned on the same line. If this is not the case, then additional parameters are 
necessary to describe the collisions from point to point. The higher the asymmetry, the 
more parameters are needed, and the problem becomes highly complicated. 

3. There is enough static friction between the rods and the floor to keep the rods from sliding 
relative to the floor. Thus the rods pivot about fixed axes. 

4. No energy is dissipated at the contact point between the rods and the floor. This condition 
is independent of the previous one; it is possible for an object to rotate about an axis and 
yet to dissipate energy. 

5. Collisions are instantaneous. This means that the time interval At during which a colli- 
sion occurs is zero. The change in total angular momentum is then AL = j^^ fdt = 0. 
That is, total angular momentum is conserved during the collision, even though gravity 
is an external force and produces a non-zero torque. 

6. The collisions are elastic. This means that energy is not dissipated during the collisions. 

7. The rods are stiff. This means that there is no deformation of the rods and thus no energy 
is converted to elastic potential energy of deformation. This condition is independent of 
the previous one; it is possible for the rods to be stiff and still dissipate energy during a 
collision. 

2.2 Definition of Symbols 

To facilitate the calculations in the next section we present our notation in advance. We label 
the rods sequentially with the numbers 1, 2, 3, . . . , /c, . . . , starting with 1 from the left end. We 
also label Ai, A2, A^, . . . , Ak, . . . , the pivot points of the rods as seen in figure [2l Then 

• 6k is the angular displacement of a rod from the vertical, and Uk = dOk/dt is the corre- 
sponding time-dependent angular velocity. 

• Vti is the initial angular velocity of the first rod immediately after it is pushed. 

• is the initial angular velocity of the k*^ rod just as it begins to move. This is of course 
the result of the collision with the {k — 1)^^ rod (for k > 1). 

• Qfk is the angular velocity of the k^^ rod just before collision with the {k + 1)*'^ rod. (The 
subscript / means 'final' or 'fallen'.) 
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• ^Ibk is the angular velocity of the k^^ rod just after collision with the (k + 1)*^ rod (b 
because the rod has just 'bounced'). 

• Pi is the angle a rod forms with the vertical at the point of collision: 

n ■ -id 

Pi = sm - . 

• Tfc is the time the /c*^ rod takes to fall from the vertical to its collision with the {k + 1)*'^ 
rod. 

In this notation, is the angular velocity of the k*^^ rod at 6^ = and ^Ifk, ^bk its angular 
velocities after it has fallen to 6* = /3i, just before and just after collision with the next rod 
respectively. 



2.3 Study of the Two-Rod Collisions 

Now examine the collision between the /c*^ and {k + 1)*'* rods. Our assumptions guarantee that 
during the collision kinetic energy and angular momentum are conserved, and that while a rod 
falls its total energy (kinetic plus potential) is conserved. These conservation laws determine 
the solution. 

Just before the collision the /c*^ rod has angular velocity Qff^ and the {k + l)*'^ rod is at rest. 
After collision the /c*^ rod has angular velocity Qhk and the {k + 1)*^ rod has angular velocity 
^Ik+i- Applying conservation of kinetic energy, 

where / = mi'^ is the moment of inertia of a rod. Therefore 

^% = ^Ik + ^Li • (2) 

Next we apply conservation of angular momentum with respect to point A^^i. Just before 
the collision the k^^ rod has angular velocity Vtjk and thus translational velocity Vk = i^jk- 
From figure 4 one can see that only the component Vk cos /3i contributes to the angular mo- 
mentum calculated around the point Ak+i, and that it does so with impact parameter £cos/3i. 
The {k + rod has no angular momentum initially. Therefore 

Linitial = m{VkCOsPi) (icOsPi) = m f Qfk cos^ Pi . 

After the collision the {k + 1)*'^ rod rotates around the point A^^i with angular velocity ^2^+1 
and thus has angular momentum IQk+i- The /c*'* rod has the new angular velocity ^2;,^, again 
around the point Ak. Therefore it will contribute an angular momentum mi'^flbk cos^/3i with 
respect to A^+i. Therefore 

L final = mf Vtbk cos^ Pi + IVtk+i . 

Conservation of angular momentum {Linuiai = L final) yields 



Vt fk cos^ Pi = Qbk cos^ Pi + Qk+i 



(3) 
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Figure 3: The series of collisions between the {k — 1) , k*^, and {k + 1) rods. In particular, 
the figure shows the state of the rods 

(a) just before the {k — if^ and k^^ rods collide. 

(b) just after the (fc — 1)*'' and k*^ rods have collided. 

(c) while the k^'^ rod rotates towards the {k + 1)*'^ rod. 

(d) just before the k^^ and {k + 1)*^ rods collide. 

(e) just after the k^^ and {k + 1)*^ rods have collided. 

(f) while the {k + 1)^^ rod rotates towards the {k + 2)*^ rod. 
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AFTER 




Figure 4: The collision between the k^^ and {k + 1) rods. 



The system of equations ([2]) and ([3]) can be solved easily for Qk+i and Qbk'- 



f- 
where 



^fc+i — /+ ^/fc , (4) 



bk 



cos^ /3i ± 1 / cos^ j3i 

Now consider the k'^'^ rod as it falls from the vertical to angle its position just before the 
collision. Conservation of total energy yields 

1 1 

-ml + mgi = ^I^}k + 1^9^ cos pi , 

Q}, = Ql+'^{l-cos(3i) . (5) 



or 



Combining equations (jlj) and ([5]) we find: 

= fM + b (6) 



where 



^9 ,2 



fiil-cos/3i] 
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Equation ([6]) is a mixed progression (i.e., a combination of an arithmetic and a geometric 
progression) and can be solved by well-known techniques (see appendix |A]) . The result is 

^ _ r2{k-l) 

Recall that Qi is the initial angular velocity of the first rod caused by the initial external push. 
We show now that /+ < 1. Since /3i 7^ 0, 7r/2, x = cos^ jSi ^ 1, 0. Then {x-ljxf > ^ + 
- 2 > ^ + > 2. From the last inequality it follows that /+ = 2/(x^ + 1/x^) < 1. 
Since /+ < 1 it follows that lim„^_|_oo /+ = and therefore 

Thus deep into the chain we find translational invariance: the initial angular velocity imparted 
to a rod by its neighbor becomes independent of position. Notice that in this limit the initial 
push given to the first rod becomes irrelevant. 



2.4 Wave Speed 



We can obtain the limiting speed of the wave by computing the time between collisions, working 
well into the chain where this becomes independent of position. 

Apply conservation of energy for the n*^ rod as it begins moving and after it falls through 
an arbitrary angle Q. This yields: 

+ mfyff = -lu'^ + mgi cos 6 . 

Setting Un = dO/dt, we can separate t from 6 and solve for the time required for the rod to 
move from 6 = to 6 = (3i: 



dt 



de 



fiH + ¥-¥cos^ 



This integral can be expressed in terms of the complete elliptic integral of the first kind K[k) 
(see Appendix [B|) : 



where = + c = ^ and /c„, = 



K{K) 



- k 



2c 
an+c' 



In the limit of large n, the time T„ approaches a limiting value 



T 



K (k) - F 



TT-/3-_ 



Kk 



where a = fi^ + ^, c = k 
V = d/T given by 



The wave therefore approaches a limiting speed 



\/a + c 



2K{k)-F{^,k) 
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A little algebra lets us write the wave speed in the scaling form of equation ([T]): 
with 



and 



t i k[K{k)-F{^,k)] 



(l-cos/3i)/2+ 2(1-/2) ■ 

Since /+ and thus k depend only on Pi = sin"^ (d/i), this G is indeed a function only of d/i, 
as required by scaling. The scaling function G is plotted in figure [51 




Figure 5: The scaling function G{d/t). 



No simple closed expression exists for complete elliptic integrals, but some insight into our 
solution comes from looking at the limit of very closely spaced rods {d <C t). In this limit 
Pi ^ d/i and ~ 1 — Pf. Using these in the above expression yields k"^ ~ 4/?f -C 1. Then 



K(k)-F{ ^^^,k 



dt 



-ft 2 



From this we find 



dli 



Thus the wave in very closely spaced rods moves very fast. 
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In a similar way one can examine the other extreme geometrical limit, d/i slightly smaller 
than unity. Put (3i = 7i/2 — \/2e. Then 

- = sin Pi = cos V2e ^ 1 — e 

while /+ ~ 4e and k"^ ^ 1 — 16e^. For k very near unity the complete elliptic integral is 
approximately K{k) ^ ln(4//c') where k'"^ = 1 — k"^ [7j. Here this gives K{k) ^ ln(l/e), which 
diverges as e approaches zero. However in this limit F is finite: 

. ~ f j?L=,„(V2+i). 

' -k^ sin2 1 Jo cost 



In the limit d ^ £ using the above approximations yields 



In 



il+V2)i£-d) 



OT V ^ y/glG{d/^), where 



g6 



i ln(l + V2) + ln(l 



as d/l increases to 1. Thus as d gets very close to £ the wave speed drops to zero. Physically 
this occurs because one rod gives a very small push to the next in line, which as a result takes 
a very long time to fall. The reader might wish to conduct a quick experiment to verify this 
conclusion or a more careful one to compare our theoretical results with experiment. Of course, 
the reader has noticed the consistency of our results with that of McLachlan et al. ([1]). 



3 Conclusions 

In this paper we have presented a set of assumptions for the propagation of the domino wave 
and we have computed the corresponding limiting speed. For simplicity we have presented the 
solution for a simplified geometry. However, the reader can easily transfer the solution to the 
case of dominoes with the shape of parallelepipeds — with appropriate adjustments of course. 
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Note Added in Proof 

After this work was completed a paper [9] with 
Cornell archives. 



somewhat similar analysis appeared on the 
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A Mixed Progression 

Consider a sequence a^, k = 1,2, . . . , with the recurrence relation 

Ofc = rak-i + b . 

This is known as a mixed progression. We want to express in terms of ai, r and b. Multiply 
both sides by r"^'^ and sum from k = 2 to n (with n > 2): 

n n 
fc=2 k=2 

In the first term on the right-hand side replace k hj k' = k — 1 and in the second replace k by 
k' = n — k. This yields 

n n—1 n— 2 

r^'^'flfc = ^ r'^~'''ak' + b^ r''' . 

k=2 k'=l fc'=0 

The first two sums have nearly all terms in common (all but the n*^ on the left and the first 
on the right). Cancelling the terms in common and evaluating the third sum yields the desired 
solution: 

1 _ ^n-l 

rj — 1 7 

an = r «! + . 

1 — r 

B The Elliptic Integral of First Kind 

The elliptic integral of the first kind [71 [8] is defined by 

f^^ dd 
nM = / , ,2-2. ' 0<A;<1. 
-'0 y\ — k'^ sm (p 

When 00 = 7r/2 this is called the complete elliptic integral of the first kind, denoted by K{k): 

'•7r/2 

K{k) 

The integral 



v^l — /c^ sin^ I 



1%) = / , , < c < a , 

Jo ya — ccosU 

can be expressed in terms of the elliptic integral of the first kind as follows. First make the 
change of variable 6 = n — 2t: 



H0o) = 2 [ 



/2-6»o/2 a/a + CCOs(2t) 
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Using the identity cos 2t = 1 — 2 sin^ t this becomes: 



7r/2-6io/2 \/{a + c) — 2c sin^ t 
2 r'^ dt 



Finally we set 



a + c 



and we rewrite the above result in the form 

dt /-/2-e„/2 



Va + c y 7o x/l - A;2 sin^ t - P sin^t 
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